Spin Waves and magnetic exchange interactions in insulating Rb .89Fei 58 Se 2 
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The discovery of alkaline iron selenide 
AFei 6+a;Se2 (^4 = K, Rb, Cs) superconduc- 
tors [1 — 5j has generated considerable excitement 
in the condensed matter physics community 
because superconductivity in these materi- 
als may have a different origin from the sign 
reversed s-wave electron pairing mechanism 
[6]-[8], a leading candidate proposed for all other 
Fe-based superconductors [9, 10J. Although 
AFei.6+ x Se2 are isostructural with the metallic 
antiferromagnetic (AF) iron pnictides such as 
(Ba,Ca,Sr)Fe 2 As 2 I12j . they are insulators 

near x — [3 — 5j and form a y/E x y/E blocked 
AF structure (Fig. la) completely different 
from the iron pnictides [3l 131 [M COB Q2J- If 
magnetism is responsible for superconductivity 
of all iron-based materials [10J, it is important 
to determine their common magnetic features. 
Here we use neutron scattering to map out spin 
waves in the AF insulating Rbo.89Fei.5gSe2. We 
find that although Rbo.sgFe! 58 Se 2 has a Neel 
temperature (Tjy = 475 K) much higher than 
that of the iron pnictides (Tjv < 220 K), spin 
waves for both classes of materials have similar 
zone boundary energies [18 20j. A comparison 
of the fitted effective exchange couplings using 
a local moment Heisenberg Hamiltonian in 
Rb . 89 Fei. 58 Se 2 , (Ba,Ca,Sr)Fe 2 As 2 |T8ti20]. and 
iron chalcogenide Fei.osTe [21 J reveals that their 
next nearest neighbor (NNN) exchange couplings 
are similar. Therefore, superconductivity in all 
Fe-based materials may have a common mag- 
netic origin that is intimately associated with 
the NNN magnetic exchange interactions, even 
though they have metallic or insulating ground 
states, different AF orders and electronic band 
structures. 

Soon after the discovery of superconductivity in iron 
pnictides |22j . calculations and experiments have found 
that electronic band structures of these materials are 
composed of hole and electron Fermi pockets near T(0, 0) 



and M(1,0)/M(0, 1) points, respectively [TU]. As a con- 
sequence, sign reversed quasiparticle excitations between 
the hole and electron pockets can induce s -wave su- 
perconductivity, giving rise to a neutron spin resonance 
at the in-plane wave vector Q — (1,0) (Fig. lc) [2"3T - 
25 . If sign reversed electron-hole pocket excitations be- 
tween r(0,0) and M(1,0)/M(0, 1) points are necessary 
for superconductivity, superconductivity in alkaline iron 
selenides should have a different microscopic origin since 
angle resolved photoemission experiments measurements 
on these materials reveal only electron Fermi surfaces at 
M(l, 0)/(0, 1) points and no hole Fermi pockets at T(0, 0) 
point 6 8. On the other hand, if AF spin excitations are 
responsible for superconductivity in Fe-based supercon- 
ductors jTUJ [25] , one would expect that spin waves in the 
parent compounds of different classes of Fe-based super- 
conductors have a common magnetic origin associated 
with superconductivity. Previous work on spin waves 
of (Ba,Ca,Sr)Fe 2 As 2 P3H2Q] and Fei. 05 Te [21] suggests 
that the NNN exchange couplings in these materials are 
similar. Since the insulating AFei.6+zSe 2 has completely 
different magnetic structure and static ordered moment 
(Fig. 1) from those of (Ba,Ca,Sr)Fe 2 As 2 and Fei.osTe 
|11) , it is important to determine if its effective magnetic 
exchange couplings are similar to these materials. 

Here we report inelastic neutron scattering studies of 
spin waves in the insulating Rbo. 8 9Fei.5 8 Se 2 with Tjv = 
475 K. Our neutron diffraction measurements on the 
sample confirmed the previously proposed Fe4 block AF 
checkerboard structure (Fig. la) [3]. Since the ferro- 
magnetic (FM) Fe4 block in the y/E x y/E superlattice 
unit cell can have either left or right chirality (Figs, la 
and lb), one expects to observe four AF Bragg peaks 
stemming from each of the chiralities. Figure lc shows 
the expected AF peaks from the left chirality in re- 
ciprocal space using the orthorhombic unit cell similar 
to that of iron pnictides [T8T - f2"0] . where they occur at 
(H a , K , L ) = (0.2 + m, 0.6 + n, L Q ); (-0.2 + m, -0.6 + 
n;L Q ); (0.6 + m, -0.2 + n,L o );(-0.6 + m, 0.2 + n, L a 
(to, n — ±2, ±4, • • •, and L Q = ±1, ±3, • • •). Considering 
both chiralities for the AF order, there are eight Bragg 
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peaks at wave vectors (H , K Q , L a ) = (±0.2 + m, ±0.6 + 
n, L a ) and (H a , K Q . L a ) = (±0.6 + m, ±0.2 + n, L a ) from 
the block AF checkerboard structure (Fig. Id), where 
the odd values of L indicate AF coupling along the c- 
axis direction [31 [THl [T7] . Therefore, acoustic spin waves 
in the AF ordered phase of Rbo.sgFex.ssSes should stem 
from these eight Bragg peaks. 

Before mapping out the wave vector dependence of spin 
waves in Rbo.89Fei.5sSe2, we first determine their overall 
energy bandwidth and the effective c-axis coupling. Fig- 
ures le and If show the background subtracted scattering 
projected in the wave vector (Q = [—1.5, K \) and energy 
plane. One can see three clear plumes of scattering aris- 
ing from the in-plane AF zone centers Q — (0, —2), (0, 0), 
and (0,2) rlu. With increasing energy, spin waves are 
gapped at energies between 75 and 95 meV (Fig. If) and 
between 150 and 170 meV (Fig. le). The zone bound- 
ary spin wave energies are around 220 meV (Fig. le). 
Therefore, in spite of the large differences in Neel tem- 
peratures of Rbo.7eFex.eSe2 (T N = 475 K) [3j QU E], 
(Ba,Ca,Sr)Fe 2 As 2 (T N < 220 K) [TSH2D], and Fe^Te 
(Tjy ~ 70 K) [5T], their zone boundary spin wave en- 
ergies are rather similar. To estimate the AF coupling 
strength along the c-axis, we show in Fig. lg spin waves 
projected in the wave vector Q = [0.6, 0.2, L a ] and energy 
space. One can see clear dispersive spin waves stemming 
from AF positions L a — 1, 3, 5 that reach the zone bound- 
ary energy near 30 meV. 

To see the evolution of spin waves with increasing en- 
ergy, we show in Fig. 2 the two-dimensional constant- 
energy (E) images of spin waves in the [H ,K ] plane 
for various incident beam energies (Ei). From their c- 
axis dispersion (Fig. lg), we know that spin waves in 
Rbo.89Fei.58Se2 are three-dimensional similar to that in 
(Ba,Ca,Sr)Fe2As2 fTMSO] and center at AF wave vectors 
Qaf = {H , K„, L ) = (±0.2 + to, ±0.6 + n, L)/(±0.6 + 
m, ±0.2 + n,L) with L Q = ±1,±3, •■■ rlu. For an en- 
ergy transfer of E = 10 ± 2 meV (above the anisotropy 
gap of E = 8 meV, see supplementary information), spin 
waves are peaked at the expected eight AF Bragg posi- 
tions Qaf around Q — (0, 0, ±1) rlu as shown in Fig. 2a. 
Upon increasing energies to E = 26 ± 2 (Fig. 2b) and 
30 ± 2 meV (Fig. 2c), spin waves from the two chiralities 
centered around the Qaf positions become apparent and 
increase in size with increasing energy. The two spin wave 
rings from the left and right AF chiralities (Figs. la-Id) 
meet near E = 45 ± 3 meV (Fig. 2d). At E = 55 ± 3 
meV, the overlapping spin waves from both AF chirali- 
ties still form rings around the Qaf positions (Fig. 2e). 
Spin waves have evolved into broad rings centered around 
(if , K 0l L ) = (±77i, ±n, L ) at E — 70±3 meV as shown 
in Fig. 2f, just before disappearing into the 75 < E < 95 
meV spin gap (Fig. If). Upon re-emerging from the spin 
gap at an energy transfer of 110 ± 10 meV, spin waves 
form transversely elongated ellipses centered at the wave 
vectors Q = (±1, 0)/(0, ±1) (Fig. 2g), identical to the 
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FIG. 1: The AF spin structure, reciprocal space, and c- 
axis spin waves of the insulating Rbo.89Fei.5gSe2- Our neu- 
tron scattering experiments were carried out on the ARCS 
chopper spectrometer at the Spallation Neutron Source, Oak 
Ridge National Laboratory. We co-aligned 2.7 g of single 
crystals grown by self-flux (with mosaic of ~6°). The in- 
cident beam energies were Ei — 80, 140, 250, 440 meV, and 
mostly with Ei parallel to the c-axis. Spin-wave intensities 
were normalized to absolute units using a vanadium stan- 
dard (with 30% error). We define the wave vector Q at 
(qx,q y ,q z ) as (H ;K ;L ) = (q x a /2Tr;q y a /2TY;q z c /2n) rlu, 
where a = 5.65 and c = 14.46 A are the orthorhombic cell 
lattice parameters. The AF spin structures are shown for (a) 
left and (b) right chirality. The v§ x v§ superlattice struc- 
ture is marked as grey with lattice parameter a s = 8.933 A. 
The orthorhombic lattice cell is shaded green. The effective 
NN, NNN, NNNN exchange couplings are marked as Ji/J[, 
J2/J2, and J3/J3, respectively, (c) The [H ,K ] reciprocal 
space with the expected AF Bragg peaks from the left chiral- 
ity. The green squares show nuclear Bragg peak positions, (d) 
Expected Bragg peaks for both chiralities. (e,f) Spin waves 
projected onto the K -E plane with H integration from -2 
to -1. The scattering were measured with Ei = 440, 250 meV, 
respectively, (g) c-axis spin-wave dispersion projected on the 
L-E plane with H integration from 0.5 to 0.7 and K inte- 
gration from to 0.4. The solid line is the calculated c-axis 
dispersion using effective exchange couplings discussed in the 
main text. 



AF ordering wave vector of (Ba,Ca,Sr)Fe2As2 
Finally, at E = 200 ± 20 meV, an energy well above the 
150 < E < 170 meV spin gap, spin waves move into wave 
vectors Q = (±1,±1) (Fig. 2h), almost identical to the 
zone boundary spin waves for BaFe2As2 [IS] and Fei.osTe 
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FIG. 2: Wave-vector dependence of spin-wave excitations 
at different energies in the [H ,K ] scattering plane for 
Rbo.89Fei.5gSe2 at 10 K. Spin wave excitations at energies (a) 
E = 10 ± 2; (b) E = 26 ± 2; (c) E = 30 ± 2; (d) £ = 45 ± 3; 
(e) B = 55 ± 3; (f ) E = 70 ± 3; (g) £ = 110 ± 10; and 
(h) E = 200 ± 20 meV. (a)-(c), (d)-(f),(g),(h) were obtained 
with Ei = 80, 140, 250, and 440 meV, respectively, along 
the c-axis. The vertical color bars indicate intensity scale in 
mbarns /sr / me V / f . u . 
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FIG. 3: Spin-wave dispersions of Rbo.89Fei.5sSe2 and fits 
using Heisenberg Hamiltonian. Spin-wave dispersions ob- 
tained by cutting along high-symmetry directions marked in 
the right panels for (a) highest energy optical energy band; 
(b) medium energy optical energy band; and (c) acoustic spin 
wave mode. The blue solid lines show fits with J3 > 0, while 
the pink solid lines are fits with J3 = 0. (d) The energy de- 
pendence of the local susceptibility and our model calculation 
of the local susceptibility. 



We use a local moment Heisenberg Hamiltonian with 
the effective nearest (NN or J±, J[), next nearest (NNN 
or J 2 , J 2), and next next nearest neighbor (NNNN or 
J3,Jg) magnetic exchange couplings (Fig. la) to fit the 
observed spin- wave spectra [TJ |23 HEl OH IS]- To ac- 
count for the ~8 meV low-energy spin gap, we add a 
spin anisotropy term J s to align spins along the c-axis 
(see supplementary information). There are 8 spins in 
each magnetic unit cell (Figs, la and lb), therefore we 
should have four spin wave bands in the Brillouin zone. 
From Figs. 1 and 2, we see that spin waves exist in three 
separate energy ranges: the lowest branch starts from 
~9 meV to ~70 meV, second from ~80 meV to ^140 
meV, and the third branch from ^180 meV to ^230 meV. 
The high quality of the spin- wave data allows us to place 
quantitative constraints on effective exchange couplings 
in the Heisenberg Hamiltonian (see supplementary infor- 
mation). While the low-energy spin waves between ^9 
meV to ~70 meV are acoustic mode arising mostly from 
AF interactions of the FM blocked spins, the two other 
branches of excitations are optical spin waves associated 



with exchange interactions of iron spins within the FM 
blocks [1, 28, 3TJH3I]. We have attempted, but failed, to 
fit the entire spin wave spectra using only the effective 
NN and NNN exchange coupling Heisenberg Hamiltonian 
(see Fig. 3 and supplementary information). For spin- 
wave fits that include the NNNN exchange coupling J 3 , 
we find that the low energy spin wave band (acoustic 
band) depends mainly on J'^J'^ J3, and J c (the effective 
c-axis exchange coupling), but not J\ and J2. The sec- 
ond band depends on the J2 heavily and the top band is 
mainly determined by J\. 

For simplicity, we consider each FM block with 4 
aligned spins as a net spin S e tt. They interact with each 
other antiferromagnetically (via Jeff) to form a cuprates- 
like AF spin structure. There is one spin-wave band for 
this effective block-spin Heisenberg model, which has an 
analytical form for spin- wave dispersion (see supplemen- 
tary information). By comparing the J e ff Heisenberg 
Hamiltonian with those of the J\- J[- J2- J' 2 - J3- J3 model, 
we find that spin waves in the first band can be approx- 
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FIG. 4: Calculated wave- vector dependence of the spin waves 
in the [H , K ] scattering for identical energies as that of Fig. 
2. The instrumental resolution is convoluted with the Heisen- 
berg Hamiltonian. 



imately described by the </ e // Heisenberg Hamiltonian, 
where J e ffS e ff = (A + 2 4 + 2J 3 )S/4 is -17 meV. This 
suggests that the low energy band is mainly determined 
by J(,J 2 , and J c . Physically, the lowest energy band 
corresponds to the block spin waves where the 4 spins 
fluctuate in phase and resemble a single spin. Only at 
high energies, the relative motions within the blocks can 
be excited, which correspond to the two high energy op- 
tical modes. Thus the high energy bands are basically 
determined by the intra-block couplings J\ and Ji- 

To quantitatively determine the spin-wave dispersion, 
we determined the measured dispersion from a series 
of high symmetry scans through the (H , H a , L a ) and 
(H . 1/2H D — 1/2, L) directions, where L D was integrated 
to improve counting statistics. Figures 3a-3c summa- 
rize the dispersion of spin waves along the marked di- 
rections on the right panels. For the low-energy acous- 
tic mode, we find a spin anisotropy gap below 8 meV 
and counter propagating spin waves for energies above 
30 meV (Fig. 3c). The two high-energy optical spin- 
wave modes are essentially dispersionless. The blue and 
pink solid lines show Heisenberg Hamiltonian fits to the 
dispersion curves with and without J 3 . The final fitted 
effective magnetic exchange couplings for spin-wave dis- 
persions are S J 1 = -36 ± 2, SJ[ = 15 ± 8, SJ 2 = 12 ± 2, 



SJ' 2 = 16 ± 5, SJ 3 = 9 ± 5, J 3 = 0, SJ C = 1.4 ± 0.2, and 
SJ S = 0.44±0.1 meV (see supplementary information for 
fits with other parameters). Figure 3d shows energy de- 
pendence of the observed local susceptibility [I] and our 
calculation using the fitted parameters. We see that the 
calculated local susceptibility agrees quite well with the 
data. To further compare the data in Fig. 2 with calcu- 
lated spin waves using fitted effective exchange couplings, 
we show in Figure 4 the two-dimensional spin-wave pro- 
jections in the [H Q , K ] plane convoluted with instrumen- 
tal resolution. The calculated spin-wave spectra capture 
all essential features in the data. 

For a Heisenberg model with spin S, the total mo- 
ment sum rule stipulates M — (gfJ,B) 2 S(S + 1). For 
irons in the 3d 6 electronic state, the maximum possible 
moment is gS = 4 ^s/Fe for g = 2, giving M = 24 
/i^/Fe. Based on absolute spin wave intensity measure- 
ments in Fig. 3d, the sum of the fluctuating moments 
below ~250 meV is (m 2 ) w 16 ± 3 (J? B /Fe. If we assume 
that the ordered moment is on the order of ^3 /is/Fe 
[31 [TBI HZ], we see that the total moment sum rule is 
exhausted for magnetic scattering at energies below 250 
meV. Therefore, spin waves in insulating Rbo.76Fei.63Se2 
can be regarded as a classic local moment system where 
a Heisenberg Hamiltonian is an appropriate description 
of spin-wave spectra. 

It is instructive to compare the effective magnetic ex- 
change couplings in different iron-based superconductors. 
First, comparing Rbo.89Fei.5sSe2 with Fei.osTe [3JJ, we 
note that although their static AF orders have com- 
pletely different structures, these two iron chalcogenides 
are very similar in terms of the values of their effective 
exchange couplings. Both of them have: (i) large FM Ji 
(or Jia), (ii) large anisotropy between the two NN cou- 
plings Ji(Jia) and J[ (or Jib), (in) AF NNN couplings 
and small anisotropy between two NNN couplings J2(or, 
J 2a ) and j r 2 (or J 2 fc), and (iv) significant AF NNNN cou- 
plings J 3 . Therefore, the presence of the iron vacancy 
ordering in Rbo.8gFei.5sSe2 reduces magnetic frustration 
and stabilizes the blocked AF structure, but does not 
change the local magnetic exchange couplings strengths 
as compared to Fei.osTe [21]. Second, comparing iron- 
chalcogenides to iron-pnictides, we find that there are im- 
portant differences as well as essential common features: 
the differences include the large FM J\ a and significant 
AF J3 in iron-chalcogenides against the large AF J\ a and 
negligible J3 in iron-pnictides, respectively, and the com- 
mon features include the large anisotropy of NN exchange 
couplings and similar AF NNN couplings. While the NN 
exchange couplings vary significantly according to the 
spin configurations between the corresponding two NN 
sites in the magnetically ordered states, the AF NNN ex- 
change coupling remains almost uniform amongst differ- 
ent classes of materials even though their AF structures 
can be quite different. This is consistent with the idea 
that J2 is mainly determined by a local superexchange 
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mechanism mediated by As or Se/Te [33 . Regarding the 
microscopic origin of superconductivity, the difference 
between the NN exchange couplings of the two classes 
of materials suggests that the NN FM exchange coupling 
cannot be responsible for superconductivity since elec- 
tron pairing is in the spin singlet channel [31], which is 
not allowed by the FM coupling. However, the similarity 
on J 2 in both classes of materials suggests that if su- 
perconductivity in all Fe-based materials has a common 
magnetic origin, it must be intimately associated with the 
NNN magnetic exchange interactions, likely resulting in 
a s-wave pairing symmetry [3SJ. 
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SUPPLEMENTARY DATA 



In addition to the spin wave data presented in the main text, we have taken triple- axis spectrometer measurements 
on HB-1 at Oak Ridge National Laboratory to determine the low-energy spin anisotropy gap. Before showing the 
results, we note that although the scattering cross section is related to the dynamic structure factor S(Q,E), it is 
proportional to the imaginary part of the dynamic susceptibility \"{Q,ui) if the temperature is much lower than the 
lowest energy spin waves. Theoretically one has S(Q,E) = 1/(1 - cx.p(- E / (k B T)))x" (Q , E). If ksT « E as is 
the case of the experiment, one has S(Q,E) oc x"{Q^E). Figure 5(a) shows x"(Q,E) at Qaf = (0.6,0.2,3), which 
clearly establishes the anisotropy spin gap of ^8 meV. Constant energy scans at 5 meV and 10 meV shown in Fig. 
5(b) confirm the presence of the spin gap below 8 meV. To further demonstrate the presence of spin gaps around 80 
meV and 160 meV, we show in Figs. 5(c)-(e) constant energy cuts for energies of E = 74 ± 4, 82 ± 4, and 90 ± 4 meV, 
respectively. There are clearly no magnetic scattering near E = 82 ±4 meV [Fig. 5(d)]. Figures 5(f)-(h) show similar 
constant-energy images at E = 140 ± 10, 155 ± 15, and 195 ± 15 meV. The scattering near E = 155 ± 15 meV are 
featureless, confirming the presence of a spin gap at this energy. 



MODEL HEISENBERG HAMILTONIAN 



The model we use to understand the magnetic excitation is a quantum spin model with up to third nearest neighbor 
(NNNN) exchange in the a&-plane, nearest neighbor (NN) exchange along the c-axis and a single ion anisotropy term, 
i.e., 

H = H ab + H c + H s , (1) 

where 

H c = J c S S r ■ S r+2 ,, (2) 



J. 



2 ) 



and H a b is given in Ref. pQ. To solve the Hamiltonian, one can use standard linear spin wave approach. A generic 
position of the spin is given by 

r = mli + nl 2 + d i; (3) 

where m, n are integers and 

li = (2x - y)/VE, (4) 
la = (x + 2y)/75, 
di = 0, d 2 = x, d 3 = x + y, d 4 = y. 

The Holstein-Primakoff transform (truncated) of the spin operators is given by 
For m + n =even: 

S+{r) = y/2Sai(R), (5) 
S-(t) = V2Sa\{K), 
S z (r) = S-ot(R)oi(R); 

For m + n =odd: 

S+(r) = 725al(R), (6) 

S_(r) - V2S ai (IL), 

S z (r) = -5 + at(R)oi(R). 

Define ip^(k) = (a\(k),a\(k),a 3 (k),a 4 (k),ai(—k),a2(—k),a3(—k),a,4 : (—k)), and we have 
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FIG. 5: Triple-axis spectrometer data obtained on HB-1 and additional spin-wave images near spin gaps on Rbo.sgFei.sgSea. (a) 
Constant-Q scan at the AF wave vector Qaf = (0.6, 0.2, 3) rlu with background subtracted, and corrected for Bose population 
factor. There is a clear spin gap below E = 8 meV. (b) Constant-energy scans across the AF wave vector at E = 5 meV 
and E = 10 meV. The data confirm the presence of a spin gap at 5 meV. Spin wave images in the (H , K ) plane for energy 
transfers of (c) E = 74 ± 4; (d) 82 ± 4;(e) 90 ± 4; (f) 140 ± 10; (g) 155 ± 15; 195 ± 15 meV. There are clearly no spin wave 
excitations at E = 82 ± 4 and 155 ± 15 meV. 



A(k) and B{k) are four-by-four matrices, defined by 
A(k) = S 



( E J ie ik * 
■ E Q 



J' z e~ i2k " 



J l€ lky 

Eq 



J. 



J x e lk 

-ik x -\-ik 



\ 



+ 4e 
Jie jfcx 
E() 



-2ik v 



(8) 



/ 2J c cos(fc z ) J> 2 e- ik *+ ik v + J 3 e 



B(k) = S 



j/g iky 



V 



2 J c cos{k z ) j^g-ikx-iky + j 3e 2ik x 

2J c cos(k z ) 



J<2& ffi iky 



J[e ik 



+ J 3 e 2lk * \ 



J' 2 e lK 



(9) 



2J c cos(fc z ) J 



where E = — (2Ji + J% — J[ — 2Jj — 2J 3 + J 3 — 2 J c — J S )S. The lower triangle elements are suppressed because both 
matrices are hermitian. 

We use equations of motion to solve this Hamiltonian. 



dip{k)/dt 



A(k) B(k) 
-B(k) -A(k) 



(10) 



Solving this eigenvalue problem for each k, we have 



H= E (7}(*)7i(fc) + l/2)wi(fc), 

i=l,2,3,4;fc 



(11) 
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FIG. 6: (a) Energy dependence of imaginary part of local susceptibility for the three different exchange parameter sets. (b,c) 
Dispersion curves for the three different exchange parameter sets as discussed in the text. 



and 



Oi(fc) = Y,U lJ {k) l3 {k) + Vy(fc)7}(-*)- 



(12) 



The differential cross section of inelastic neutron scattering can be expressed in terms of the spin wave dispersion and 
wave functions: 



^f^i s) — + "B^i T)) | Ui a (q) + V* a (—q)\ 2 D(uj, 0J a ). 



(13) 



In the above expression, I (u},q) includes all factors of experimental resolution extracted from information of each 
detector, fie(w, T) is the Bose factor and D{uj,uj a ) is the harmonic oscillator damping given by 



7T (oj 2 -lu 2 ) 2 + 4T(uj) 2 uj 2 ' 



(14) 



The damping strength T(uj) is approximated by a linear function of energy whose explicit form is to be fitted. Our 
fitting is based on so far the most general spin model with all symmetry allowed exchanges up to NNNN. A failure 
of this model in understanding the data would mean that the observed excitations cannot be explained by a local 
moment picture and the effect of itinerant electrons must be seriously considered. 

FITTING CONSTRAINTS 

The high quality of the data allows one to place quantitative constraints on parameters in the model. The data 
shows that the excitations exist in three separate energy ranges. The lowest branch starts from ~ 9 meV to ~ 70 
meV, second from ~ 100 meV to ~ 140 meV and the third branch from ~180 meV to ^230 meV. The low energy 
part of the first branch can be fitted very well by the form 



e(k) 



£ k 2 



(15) 



10 



(a) Experiment (b) J set (1) 



(c) J set (2) (d) J set (3) 




-2-1012 -2-1012 -2-1012 -2-10123 
H (r.l.u) H (r.l.u) H (r.l.u) H (r.l.u) 

— Jset(1) J set (2) J set (3) 



12 



"I 1 1 ' 1 — 

E = 97.5 ±7.5 meV 



T ' 1 ' 1 ' 

E = 9 7.5+7.5 meV " 




T 1 1 ' 1 — 

sE = 97.5+7.5 meV 



0.0 0.5 1.0 0.0 0.5 1.0 0.0 0.5 1.0 

[-1, K ] (r.l.u) [-1, K ] (r l.u) [ -1, KJ (r.l.u) 



6 12r 




-2.0 -1.5 -1.0 -0.5 -2.0 -1.5 -1.0 -0.5 -2.0 -1.5 -1.0 -0.5 
[H o ,0] (r.l.u) [H o ,0] (r.l.u) [ H , ] (r.l.u) 




0.0 0.5 1.0 0.0 0.5 1.0 0.0 0.5 1.0 

[-1,K ] (r.l.u) [-1, Kq] (r.l.u) [-1,K ] (r.l.u) 

FIG. 7: Calculated dynamic structure factor and their comparison with Heisenberg Hamiltonian with different exchange pa- 
rameters, (a) Constant energy cut of data at E = 102.5 ± 7/5 meV projected onto the (H ,K ) plane. (b,c,d) Calculated 
dynamic structure factor S(q,w = 102.5 ±7/5) projected onto the (H ,K ) plane for three different exchange coupling parame- 
ters. (e-m)Cuts along different directions and their comparison with spin wave calculations in three different exchange coupling 
parameters. 



with Vg Xp = 300 meV • A and A e s xp = 9 meV. At the propagation vector of the ground state Q = (0.6, 0.2, 1) rlu (in 
the orthorhombic basis), energy has k z dispersion, and the band top is about E% xp ~30 meV. All these values have 
analytical expressions in the spin wave model. The anisotropy gap (bottom of the first branch) is 



A s = S^J S {2J[ + AJ> 2 + 4J 3 + 4J C + J s ). (16) 
The top of the first band is reached at Q a = (0.2, 0.4, 0) rlu with 

E lt = 2S[2jf + J[{J' 2 + J 3 -J^ + J c ) + (J 2 + J 3 - 4)(4 + J 3 -J' 3 + 2J C ) - J X (JJ + 2(J' 2 + J 3 -J' 3 + J c )) (17) 



*4Jf + J? J* + U\{,V 2 + J 3 ~J 3 + J C ){J[ + J! 2 + J 3 -J' 3 + Jc) - 4Jj 3 (J( + 2{J' 2 + J 3 -J 3 + J c ))]5. 
Without single ion anisotropy, i.e., J s = 0, the spin wave velocity is given by 

v s - \[\s{[MJ[ + 2{J' 2 + J 3 - J£) + J[(J 2 -J' 2 -J 3 + J' 3 ) + 2{J 2 {J' 2 + J 3 - J' 3 ) + 2J 3 (J' 3 - J' 2 ))] (18) 

(J{ + 2{J' 2 + J 3 + J c ))/(Ji -J[ + J 2 -J 2 -J 3 + J' 3 )} 1/2 . 
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The expression with J s ^ is also available but too lengthy to be placed here, and interested readers can request it from 
the authors. The second branch actually contains two close spin wave bands. The branch starts at Q = (0.3,0.1, 1) 
rlu with energy E 2 b, whose expression is again too lengthy to be published. The second branch ends at T = (0,0,0) 
point with 

E 2t = S^J (2Ji - 2 1[ + 2 J 2 - 2J' 2 - 2 J 3 + 2J' 3 - J s )(2Jx + 2 J 2 - 2J' 2 - 2J 3 + 2J' 3 - 4 J c - J s ). (19) 
The highest branch starts at T point with 

E 3b = (4Ji - 44 - 4J 3 - J s )(4,h - 2J[ - 4J C - J s ), (20) 

and ends at (0.2,0.4,0) with 

E 3t = 2S[2Jl + J[(J 2 + J 3 -J 3 + J c ) + (J 2 + J 3 - J 3 )(J 2 + J 3 -J' 3 + 2J C ) - J^Ji + 2(J 2 + J 3 -J' 3 + J c )) (21) 
+^4Jf + J[ 2 JZ + AJl{J' 2 + J 3 -J' 3 + J C )(J[ + J 2 + J 3 -J 3 + J c ) - 4J3(J( + 2(J' 2 + J 3 -4 + J c ))]i 

The band top along the c-axis is reached at (0.6,0.2,0) with 

E c = Sy/[2(J{ + 2J' 2 + 2J 3 ) - J S ](4J C + J„). (22) 

Based on the data and considering the effect of large damping at high energies, we have for the above quantities the 
following constraints: 



A, 


. - \exp 
s 


= 8-12 meV, 


v s 


= v e s xp 


= 250 - 300 meV • A, 


Eu 


= 60 ~ 


75 meV, 


E 2 b 


= 90- 


110 meV, 


E 2 t 


= 110- 


- 130 meV, 


E 3 b 


= 180- 


- 200 meV, 


E 3t 


= 200- 


- 220 meV, 


E c 


= 25- 


30 meV. 



FITTING PARAMETERS 

The above constraints give a very narrow range of parameters, we can further constraint possible exchange constants 
so that a quantitative fit to the data shown in the paper can be found. In this section we discuss what elements are 
indispensable to our fittings. 

We first emphasize that a proper fitting should have J 3 > and J[ > (antiferromagnetic) . To see this, we compare 
the following possible parameters since they can all approximately describe the data: 

(1) S,h = -36, SJ[ = 15, SJ 2 = 12, SJ' 2 = 16, SJ 3 = 9.5, SJ 3 = 0, SJ C = 1.4, SJ S = 0.44 meV. (2) S.h = -36, 
SJ[ = -5.7, SJ 2 = 13.4, SJ' 2 = 22.4, SJ 3 = 14.2, SJ' 3 = 0, SJ C = 1.4, SJ S = 0.44 meV. (3) SJx = -36, SJ[ = 10, 
SJ 2 = 11, SJ' 2 = 28.7, SJ 3 = 0, SJ' 3 = 0, SJ C = 1.4, SJ S = 0.44 meV. Figure 6 summarizes the calculated x"M and 
spin wave dispersions for all three sets of parameters. From the calculation, we see that all three parameter sets give 
similar local susceptibilities, and therefore cannot be distinguished based on x"{lo) alone. 

By comparing the calculated spin wave dispersion curves with data, we were able to separate which model is 
correct. Figure 6(b) and (c) shows the outcome for the three sets of exchange couplings for the acoustic and optical 
modes, respectively. We see that parameters of (1) and (2) fit the acoustic and optical data slightly better. Although 
the imaginary part of local susceptibility and dispersion curves for different exchange parameter sets are similar, 
their constant energy patterns at —110 meV are very different, which provides key clues to the choice among different 
exchange coupling parameters. In the energy range around 110 meV, several optical branches are mixed together. The 
combined spin wave intensity patterns depend sensitively on the exchange coupling parameters. Figure 7 compares 
directly the calculated patterns with the observation for the three set of exchange parameters. Clearly, the first set of 
parameters describes the data much better. This is what we have used to determine the effective magnetic exchange 
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FIG. 8: Spin-wave dispersions of Rbo.89Fei. 5 gSe2 and fits using Heisenberg Hamiltonian with three different exchange coupling 
parameters as discussed in the text. Spin- wave dispersions in the acoustic branch obtained by cutting along high-symmetry 
directions and model fits using three different sets of exchange coupling parameters (a)-(o) Cuts along the [—0.6, K ] direction 
by integrating H from —0.65 to —0.55. (p)-(u) Cut along [H o ,0.2] direction by integrating K from 0.15 to 0.25. 



coupling constants. This conclusion is further confirmed by comparing the calculated dispersion with the observed 
dispersion using the three sets of parameters as shown in Figs. 7, 8 and 9. 

As a remark, we note the important fact that the in-block NNN exchange J2 must be positive (antiferromagnetic) 
for all candidate sets of parameters. J2 has little effect on the first and the third branches of dispersion, but is strongly 
coupled to the middle branch. A ferromagnetic J 2 can push up the second branch for about 30%. This means the gap 
between first and second branches would be more than 40 meV, while in experiment it is clearly less than 30 meV. 



SUM RULE 

Here we discuss the total moment sum rule. For a Heisenberg model with spin S, the sum rule is formulated as 
Ref. g]: 

M = ^Y< f dk I d "S aa (K u) = M x + M y + M z = g 2 fi 2 B S(S + 1), (24) 
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FIG. 9: Spin-wave cuts of Rbo.89Fei.58Se2 and fits using Heisenberg Hamiltonian with three sets of parameters. Q-cuts for the 
highest branch. The H integration range in (a),(c) is from —1.05 to —0.95. Integration of the K range in (b),d) is from 0.8 
to 1.2. 



where g is the Lande factor. For free electrons g = 2. In Rbo.8gFei.58Se2, the maximum possible spin S = 2 is 
expected, which gives M Q = 24 fi B /Fe. 

The longitudinal part M z comes from the static moment (elastic) and the inelastic contribution. For our system, 
the static moment is about 3 /.ib/ Fe [3], which contributes 9 /i^/Fe. The inelastic part mainly comes from the 
two-magnon scattering process. The magnetization reduction can be evaluated as AS* = 0.5 from the static moment 
for S = 2. From Ref. [2], we can estimate the two-magnon spectral weig ht as AS (I + AS)g 2 (i 2 B ~ 3/x|/Fe, where 
the normalization factor has been chosen as 1. The spectral weight from the two-magnon process is only 1/3 of the 
elastic part, which is much weaker than the cuprates which has S = 1/2. In unpolarized neutron experiments, the 
two-magnon spectral weight is generally very hard to detect. We will ignore it in the following treatment. 

The transverse part M x + M y mainly comes from the one-magnon spin wave spectrum. According to Eq. (1) in 
Ref. [4], we can get the dynamic structure factor S(E) by removing the magnetic form factor. Then using Eq. (5) in 
Ref. [3], we can get the transverse part by integrating S(E) over the whole energy range. Experimentally we do not 
observe the neutron scattering signal above 250 meV, so we can choose the integration range from 8 to 250 for the 
inelastic signal only. We get the transverse part ~ 26±5 //^(/.u.) -1 , where f.u. means formula unit. Considering the 
formula of Rbo.89Fei.5sSe2, we divide it by a factor of 1.6. The transverse part M x + M y is evaluated as 16±3^i B /Fe. 

The total moment from our evaluation is 25 ± 5 [i B /Fe, which is very close to the expected total moment from the 
sum rule. Thus the Heisenberg model with S = 2 is an appropriate description for the insulating Rbo.89Fei.5gSe2 and 
the spin waves describe the spin dynamics very well. 
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